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SUMMARY 


The spanwise vortex-lattice arrangement is mathematically established by 
lattice solutions of the slender wing which are shown to be analogous to the 
chordwise vortex-lattice thin wing solution. Solutions for any N number of 
panels to infinity are obtained. With the optimum lattice for any N value the 
slender wing theory lift and induced drag and thin wing theory lift and moment 
are predicted exactly. For N"^, slender wing elliptic spanwise loading and 
thin wing cotangent chordwise loading are predicted, which proves there is 
mathematical convergence of the vortex-lattice method to the exact answer. 
Based on this A-k) and an A*^ planform spanwise lattice arrangements, an 
A-vortex-lattice spanwise system is developed for arbitrary aspect ratio. 

This A-lattice has the optimum characteristic of predicting lift accurately 
for any N value. 


INTRODUCTION ' : 


Growth of computer facilities has given the engineer a powerful tool for 
obtaining solutions to generalized problems. This is possible because with 
numerical or finite-difference methods the equations of a problem can be 
simplified readily to computer language. Vortex-lattice methods have been 
developed extensively for steady and unsteady pressure prediction and for 
planar and nonplanar configurations. Examples of some of the work and 
investigations in the vortex-lattice method and aspects and applications of 
this method are reported in references 1 through 20. A brief description 
of the typical vortex lattice is that the surfaces are divided in the spanwise 
and chordwise directions into panels which cover the surface with a lattice. 
The sides of the panels are parallel to the freestream and the chordwise panel 
boundaries follow the surface contour. The 1/4 chord line of each panel 
contains a bound or load vortex while the trailing vortices are at the sides 
of the panel. The boundary condition of no flow through the surface is ful- 
filled on every panel at one point located at the lateral center of the 3/4 
chord line of the panel. These panels are distributed in a uniform, and thus 
geometrically simplest mesh, referred to as a planform lattice. However, in 
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1966 an application of the lattice method of reference 5 to a swept wing showed 
that the chord loading in the panel bordering the wing leading edge was too t 

low, the spanwise loading near the wing tip was too high for engineering 
acceptability, and the net lift slightly too large. This was with ICO panels 
on the semispan, 10 chordwise times 10 spanwise. Investigation of a lattice 
mathematical model was made at that time of the chordwise panel distribution 
and later reported in reference 20. The results of this work showed that the 
loading at the leading edge panel needs a factor of about 1.128 which improved 
the loading value in that panel. However, the too high wing tip loading and 
lift was not explained. In early 1972 a mathematically rigorous spanwise 
vortex-lattice analysis was developed based on slender wing theory (part of 
ref. 19). This was mathematically analogous to the ea*"lier chordwise solution 
but more complicated. This spanwise lattice arrangement is characterized by 
a 1/4 - 3/4 rule which locates the trailing vortices inboard 1/4 of the 
planform panel span and the no flow through points inboard 3/4 of planform 
panel span from the planform panel outboard edges. This inboard shift of 
the lattice leads to solutions with less loading near the wing tip and less 
lift, which improves the loading in the above example. This example supports 
the observation that accuracy depends on the position of the panels in the 
lattice as well as density of panels. 

The objectives of the present study are to correlate and extend the work 
of references 19 and 20, to investigate the effect of three-dimensional 
planform on lattice arrangement, and to formulate a generalized vortex-lattice 
arrangement and method for three dimensional wings. 


SYMBOLS 


A 

aspect ratio 



swept panel aspect ratio [eq. (52)] 


b, c 

wing span, wing chord 


'L, 

lift coefficient, lift-curve slope 


COi 

induced drag coefficient 



section lift coefficient, section lift-curve slope 


a 

section pitching moment coefficient due to angle of 

attack 

®Mc 

parameter of chordwise loading [eq. (38)] 


%nc 

chordwise loading factor [eq. (41) and table 5] 


G 

spanwise loading coefficient or dimensionless circulation [eq. (1)] 

®Nn 

parameter of spanwise incremental circulation [eq. 

(7)] 
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h spanwise loading gradient factor [eq. (19) and table 2] 

Nn 

■j spanwise loading factor [eq. (22) and table 3] 

Nn 

N, Nc integer number of panels on wing semispan, and wing chord respectively 

n,m integers denoting spanwise position of vortex, and downwash point 

nc,mc integers denoting chordwise position of vortex, and downwash point 

V free stream velocity 

a angle of attack 

r circulation, also Gamma function 

n lateral coordinate per wing semispan 

Ai /2 sweep angle at 50% chord line 
^ longitudinal coordinate per wing section chord 

Subscripts: 

V vortex 

w/ downwash point 

LE leading edge 


SLENDER WING OPTIMUM VORTEX LATTICE 


Physical Similarity of Trailing Vortex Sheet Flow with 
Chordwise Thin Wing Theory Flow 


The objective of this study is to do a rigorous analytical derivation 
to determine the optimum spanwise distribution of panels analogous to the 
analysis done in reference 20 for the optimum chordwise distribution of 
panels. Optimum here defines the lattice which best duplicates exact 
solutions. A physical similarity does exist between the vorticity distribution 
of the chordwise loading with the trailing vortex sheet from a finite span 
wing. This can be seen graphically in figure 1 where thin airfoil theory 
chord-load vorticity is compared with the trailing vorticity which is predicted 
by slender wing theory (refs. 21 and 22). It has often been noted that the 
mathematics of thin wing and slender wing theories have a striking similarity. 


"om figure 1, it is noted that a similarity of vorticity is obtained 
when uhe winq tip at n = 1 correlates with the wing leading edge, and the 
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midspan point at n = Q correlates with the wing trailing edge. Since this 
correlation makes the flow fields analogous, it follows that the optimum span- 
wise panel distribution is analogous to the optimum chordwise panel distribu- 
tion given in reference 20. This is subject to the condition that the 
distribution start at the wing tip and proceed inboard. Applying this 
condition and using these distribution conditions, the spanwise panel 
distribution becomes that shown in figure 2. The determination of an optimum 
chordwise panel distribution is made by two-dimensional i zing the problem to 
planar flow and thus the spanwise extent of the panel is infinite. The optimum 
chordwise panel distribution is that which yields thin wing solutions which 
most accurately duplicate the results of exact thin airfoil theory. In the 
present work, the determination of an optimum spanwise panel distribution 
will be made by two-dimensional i zing the problem to cross-sectional flow 
and thus the chordwise extent of the panel is infinite. The optimum A->0 
spanwise panel distribution is that which yields solution which most 
accurately duplicates the results of slender wing theory. 

The objective of the present paper is to apply the analytical methodology 
of reference 20 to determine the optimum spanwise panel distribution. This 
distribution should result in an exact prediction of total lift for any 
number of spanwise panels and provide spanwise loading factors. The chordwise 
panel distribution analysis is correlated with thin airfoil theory for the 
determination of optimal accuracy. In an analogous procedure the spanwise 
panel distribution analysis here will be correlated with slender wing theory 
for the determination of optimal accuracy. 


Formulation of Spanwise Lattice Matrix and Solution 
to Infinity which Satisfy All Boundary Points 


Slender wing theory equations for additional loading (ref. 22) are 


C 


L 


a 


ttA 

2 ’ 


G„{ri) 


(1 - n2)i/2; 


dG^(n) 

dn 



c c 

where G^(n) = 

Also presented in reference 22 are solutions for flap, ailerons, and all 
spanwise loadings, which can be used to evaluate lattice accuracy when a 
problem involves these types of loadings. 


( 1 ) 


By Biot-Savart law, the downwash at nwm due to an infinite extent vortex 
at Oyn is (see fig. 2) 


AF 


^ n/b 

mn - —7 

•^(n 


vn " '^wm^ 


(2) 


with vortices located at 
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* "1 


n = ± 

vn N 

boundary condition points at 


*^wni ^ 


m - 3/4 
N 


(3) 


(4) 


with equations (3) and (4) the downwash angle at is 


t (— ! ) 

“ \n - ''wm \n * "vm 


^Nn 

3N 


(5) 


with equations (3), (4), and (5) for additional loading (am=a), equation 
(5) becomes 

N 

^ " X I ^2n - 2m + f 2n + 2m - 2^ ^Nn ; m = 1 , 2 , . . . N 
n=l 

where 


(6) 


3NAr^,„ 
9 m_ - Nn 


^Nn = 


nbVa 


= - Nag. 


*Nn 


(7) 


Generalized Inversion of g^^^^ Equations 

Equation (6) represents N unknowns, g)\|n, and N Equations. An inversion 
of equation (6) means a linear simultaneous solution of N equations. Solutions 
for N = 1, 2, . . . can be obtained readily for small N and from the resulted 
series formed of g^j^, the general solution can be determined by induction. For 
N = 1 , the solution is 


9n = 1 

For N = 2, the linear simultaneous solution of two equations is 
^ ^^11 T 8 ^ 22 ’ ^21 5 

1 1 n 4. 7 „ n 54 

' ■ ■ 2 ^21 9 ^22 ^22 ' 35 


( 8 ) 


These solutions are done for higher values of N until a sequence is formed. 
This sequence is presented in table 1. Examination of the 1st column in 
table 1 shows the sequence follows the general term of 




/ 
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1. 





I 


i 




for n = 1 : 


9n1 


3N(2N - 
(2N*l)odd' 


In the second column, ratios of 9|,^2^9 nL sequence 

N = 2 3 4 5 6 


% , 32x3 3 ^ 32x1 5 

gj^l ~ 1x7 3x9 5x11 


32x24 32x35 

7xiy 95^ 


This shows the general term as [and using gj^-j from eq. (9)] 


for n = 2: 

32(N2 - 1) 3x32(N+1)(N-1)(N)(2N-5)q^^! 

= T^N“r3lT2r:Tr " (2 n+3)^^^i 

Similarly, ratios of 9|^j3/g|^2 sequence 

N = 3 4 5 6 

M = 5^x5 52x12 52x21 52x32 

^N2 1^^^ 3x13x4 5x15x4 7x17x4 

This shows the general term as [using gj^^ 69- (10)] 

for n = 3: 

52(N2-22)gj^2 3 x32x52(N+2)(N+1)(N)(N-1)(N-2)(2N-7)^^j! 

" 22(2N-5)(2N+5) " 22(2N+5)'qJ^! 

These sequences of gj^^ , gj^j2» and g^j^ show that 

9n4 

■ 32(2N-7)(2N+7) 

^ 3x32x52x72(N+3)(N+2)(N+l)(N){N-l)(N-2)(N-3)(2N-9)^^jj 

3^x2^(2N+7)^^^! " 


(9) 


(10) 


( 11 ) 


( 12 ) 


Then in general the recurrence is [from first equalities of eqs. (10), (11), 
and (12)J 
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g„ ,n-^/2^9 N2-(n-l)2 g^, ^ 

Nn - ( n_i N2-(n-l/2)^ ^n-1 


The general term for net values of g,. is [from second equalities of 
eqs. (10), (11), and (12)] 

^ 3(N+n-l)![(2n-l) ..,]2 (2N-2n-l) ... 


(n-1 )!L(n-l )!J^ 




Equation (14) is an exact mathematical inversion of the matri . type represented 
in equation (6). The equation (14) solution is exact for a.' interger value 
of N, including N equal to infinity. 

The odd factorial (2n-l) .^i can be converted to direct 'actorials by the 
relation 

(2n-D! = (2n-l)^^^! (2n-2)^^! 

where (2n-2)g^! = (2n-2) (2n-4) (2n-6) .... = 2'’’^n-l)(n- )(n-3) 

.... = 2""'' (n-1)! 

then (2n-l)ojj! = (2n-l ) 1/2"''’ (n-1 ) ! (15) 

Using equation (15) in equation (14), then Qj. in terms of conventional 
factorials is 


q„ 3n2(2N-2n)! f(2n)!(N+n)n 2 

■ (N+n)(2N+2n)l [(n! )^(N-n)! J 


Factorials suggest Gamma functions. Extensive relations of Gamma 
functions and tables are presented in reference 23. In terms of Gamma 
functions equation (16) becomer 

9mo - 3n2 r(N-n+l/2) r(N+n-t-l) f r(n+l/2) ~l 2 


■ Tr(N+n) r(N-n+l) r(N+n+l/2) L on^rj 

The gfjn function is thus expressed in three forms given by equations 
(14), (16), and (17) respectively. 


Gradient of Spanwise Loading 

The spanwise loading gradient is [using eqs. (3) and (7)] 
^a' ^’^vn) = " I ^Nn 


k 


I 


i 



whF.re 9 i^n is given in equation (16). The spanwise position, Hyn. is given by 
equation (3). Comparable slender wing theory values of Gct'(nyn) are obtained 
by replacing n by nvn i>^ gradient function given in equation (1). Com- 
parison of equation (18) with the slender wing theory value of loading gradient 
shows that the lattice-method loading gradient requires a factor. This factor 
can be formulated accurately by this slender wing solution. 


The loading gradient factor is defined as 


h|^^ _ slender wing theory 




(19) 


Numerical values of h^n computed from equations (19) and (18) are presented 
in table 2. These factors are very near unity except at the wing tip 
region. Included are values for N-^ determined in a following section. 


/ 


t 




1 

. I 


Spanwise Loading 

With the vortex position set by equation (3) the AGfjn extends from nyn = 
(n - l/4)/N to(n - 1 - l/4)/N. The middle of this segment is at nn = 

(n - 3/4)/N, that is at the same spanwise position as the boundary condition 
points given by equation (4). The spanwise station of the loading will be 
assumed to be at the middle of the segment, that is, at nn- symmetrical 
loading, the loading at wing center is constant in the range 
-(1 - 3^)/N - n - (1 - 3/4)/N, and the middle of this segment is at = 0. 

From equation (7) the loading at the n=N segment is 


®aNN ^ ^®aNN " 3N ^NN 


At N-1 segment 
Thus at the m segment 


®aN,N-l ^^aN,N-l " 3N 


aNn 




^Nm 


m=n 




( 20 ) 


where g|\jni "'s given in equation (16) and in table 1, but with n=m numbers. 
These loadings are at the spanwise stations 


n 


n 


n - 3/4 


N 

0 


; for n > 1 
; for n = 1 


(21) I 


Comparable slender wing theory values of G(j(nn) are given by equation (1) with 
nn of equation (21). Comparison of equation (20) with the slender wing theory 
value of spanwise loading shows that the lattice-method loading requires a 
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factor. This spanwise loading factor is defined as 

iK,_ _ slender winq theory ,r, 

"" 

Numerical values of i|^n, computed from equation (22) are presented in table 3. 
These loading factors are very near to unity. Included are values for N-^. 


Lift-Curve Slope and Induced Drag 

The lift coefficient is the integration of the loading coefficient, then 

1 

C|_ = A G(n) dn (23) 

Since A6|jin is constant between nvn "’’vn* integration for lift is 

a summation of the pyramid layers of 2nyn '^^n* equation (23) 

N 

n=l 

With the AGrvNn given in equation (7) and nvn in equation (3), then equation (24) 
becomes The lift-curve slope is 

N 

‘"'-a ' Ir Z (n - l/4)g„„ (25) 

n =1 

where gwn is given in equation (16) and in table 1. With equation (16) 
inserted into (25) the lift-curve slope becomes CLc = ttA/ 2 for any value of 
N. This compares with the slender wing value given in equation (1). 

In a following section for N=“ it is proven that Cl^ = ttA/2 for 

The induced drag coefficient is given by (for constant a^) 

N 

Si = ^ _/a.G(r,)dn = I ^ VSn^n 

n=l 


where an is that in equation (5) but with changed position of n and m. Since 
an=a by the conditions of equation (6), and using equation (24), then Coi = 
C[_a/2. It was shown by equation (25) that Cl = irAa/2 for any value of N, 
then 
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which is then also valid for any value of N. This is identical to slender 
wing induced drag. Therefore, it is concluded that this lattice and boundary 
point distribution results in exact integrations for lift and induced drag for 
any interger N. This exactness was not unexpected since the mathematics 
is similar to the chordwise solution (ref. 20) in which the first harmonic 
solution (elliptic chord loading) integrates exactly. 


Solutions for N Approaching Infinity 

Detailed mathematical derivations are developed in reference 19. Here 
the results are a digest of the mathematics in that study. The primary 
purpose for exploring solutions at N'^ is to prove that a finite lattice- 
method solution does mathematically converge to the exact solution. In this 
problem the exact solution is elliptical spanwise loading evaluated from 
slender wing theory. Mathematical proof is needed that the integration for 
lift remains exact as N-^. Asymtotic values of the factors hNn 
and if^n as are useful in the tables of these factors. As t!ie 
functions become so unique and manageable the reader will find a mathematical 
exciting experience. 

Equation (16) can be used to mathematically prove that spanwise 
distribution predicted from finite element loading methods do converge on 
the exact answer as the number of elements are increased. 

Using the relation for large factorials 


N! iq: (2,)'/2 


(26) 


then equation (16) becomes (except at the point n=N, i.e. at Ti=l when N=“) 


^Nn = 


3n (N n)^^^ 
7 t(N + n) (N - n) 


W 


3^ 

’'(' - jij-2) 


Combining equations (3l, (18), and (27) results in 

-n. 


G ’(n„ = 


vn 


a ' vn 




1/2 


(27) 


(28) 


Equation (28) is the same as the loading gradient given in equation(l) by 
slender wing theory. Thus it is proven that as lattice panels are increased, 
the solution converges to the exact loading. 

With the relation tha.*^ 
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r(N tj/il n:!; Wto 

r(N + i) (N + 1/4)1/2 

then equation (17) can be expressed in the forms 


n * finite 

9n - 3n2 fr(n + 1/2) ^ 
■ ttN r n + 1) , 


N-~ 

N-n = finite 


^N,N-n 


(30) 

3Ny^ T(N - n + 1/2) 
r(N - n +1) 


Then for N'^ the spanwise loading gradient factor of equation (19) becomes 


n = finite 
^Nn = 


n - 1/4 rr(n + 1) 
~ 1 \ Lr(n + 1/2) 


N-n = finite 
^N.N-n = — 


r(N-n + 1 


(N-n + l/4)‘''^r(N-n + 1/2) 


With the use of the series summation 


I! • 2(N - m * 1/2) 


r(N - m 4 1/2) 
riN - m + ) ) 


the spanwise loading factor of equation (22) for N»®becomes 

(N - n + 3/4)^^^ r(N - n 1) 

^N,N-n ■ (N - n + 1/2) r(N - n + 1/2) 

Numerical values of equations (31) and (33) are listed in tables 2 and 3. 
Examination of these equations indicates simple relationships. These 
relationships extend to chordwise loading factors developed in a later 
section and can be expressed in one equation as follows: 


for N = N^ = “ T /o 

C h ' ^ 1 

h ^Nn _ 'N.N-n+1 

N,N-n - ,75- - 7 ^ 




(1-1/I6n2)'^ [l-l/16(N-n+l)2]'^^ [l-l/16(N,.-ti,.-l)2] 


For the zero condition 


Nc.nc-1 


’N,(N-n)=0 " \.(n^-l)=0 " 


= 1.128379 
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For N-n = 1 

^N,(N-n)=1 " \.(n^-l)=l 1.009253, 


= .954930, '|^_(N-n-tl)=l 


\ (N -n ■ -^”205 

C, ' C C 

For N = “ the expression for C|_ develops into the following series 
summation: “ 


A=0 


This factor of A/2 is tt, thus = nA/2. 


CHOROWISE VORTEX LATTICE 


Formulation of Chordwise Lattice Matrix and Solution 

to Infinity 

In reference 2c the lattice is distributed into equal length chord 
segments and the chcrd loading vortex is located ot the 1/4 points of each 
of the chord segment':; and boundary conditioT point at the 3/4 point of 
each of the chord segments. Then the load vortex chord station and the 
boundary condition chord station are respectively at 

r 1/A -1/4 

^ = Ik-3/4 . . = c / 

c ' c 

where c = x/c. 

The loading equations to be solved are given as. 


Z/ 2m -2 


i^-2n^+l = 1; ^ • ' * ^c 


where 


N,r„, Y„, '‘'pn. 


e„ „ . "c' "c 'nc 

W " TTCVa " TrVa 


( 
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The matrix inversion of equation (37) is obtained as a factorial function 


e.. <2Nc-V'>0dd! (2"c-3)»id! 

‘ (2"c-2"c>ev! 2'"‘=-'[(N^-dJ!(n,-l)!]' 


In terms of Gamma functions 


■\"c 


2n^(2Nj.-2n^+l) r(n^+l/2) r(N^-n^+1/2) 


Tr(2n^-1) r(n^+l) r(N^-n^+l) 

Numerical values of equation (39) are presented in table 4. 




(39) 


Chordwise Loading Factor 

The exact solution by thin wing theory gives the additional loading by 
the function [using defined in eq. (36)] 



= 2 J~^c ^l/2 ^ 2 ( 

irVa 71 ' ^ TT 'n^-3/4 


(40) 


The chordwise loading factor is defined by the ratio of equation (40) to 
equation (39). 

f„„ ^"cncthin ^"cthin r(W^-n^^l) 

"" ®N^n^ n^(2N^-2n^+l) (n^-3/4)^/^ r(n^+l/2) r(N^-n^+l/2) 

Numerical evaluation of equation (41) for a range of N^ are given in table 5. 


Chordwise Lift-Curve Slope and Pitching Moment 

Section wing lift is the integration of the chordwise loading and section 
moment is the integration of the product of chordwise loading and moment arm. 
In mathematical representation 

^0^ VF ^ 0^ fe 

The lattice loading is constant with c over the interval 

AC = [n -(n„-l)]/N = 1/N and the lattice load vortex is located at 
n^ c c '■* c c 

C^c = (Oj,-3/4)/N^ along the chord and the moment arm extends to this vortex. 
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Then with equation (38) equation (42) becomes 


c 


z 

a 






(n^-3/4) e 


"cV 


(43) 


The summation terms in equation (43) are listed in table 4 which when inserted 
into equation (43) show that the lift and moment for arbitrary N is the same 
as predicted by thin wing theory, that is 


s ' 2’’ 

a 


m 


\E 


2 ’ 


a.c 


'm 


\E 


= 1/4 


(44) 


Proof that c^^ is Exact at All Values of N 
Inserting equation (39) into (43) results in 

c 4 v' "c<2"c-V> r(N^-n^+l/2) 

Z (2n^-l) r(n^+l) r(N^-n^+l) 

Now r(n+l/2) = (n-1/2) r(n-l/2); r(N-n-l/2) = r(N-n-f-l/2)/(N-n-l/2) ; and 
r(n+l ) = nr(n) = n! , then 


c 


z 

a 


= 4 



r(n +1/2) r(N -n +1/2) 

n ! (N -n )! 
c ' c c' 


(46) 


The Gamma functions and factorials in equation (46) show ttjt this product 
term is a symmetric function which is factored by an antisymmetric term 
(Nri\)/I^. Then some of the high Oq terms cancel the low nc terms. By expanding 
the summation a new summation can be formed given by 


c 


z 


a 




c 2 



Nc 

^(^-n^) r(n^+l/2) r(N^-n^+l/2) 


(47) 
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where 


- n^) = 1 for = ^ .= 2 for ^ 

The summation term in equation (47) is the Legendre polynomial of the first 
kind Pn (cos 0)q=o (see p. 36, ref. 24) and Pf^(l) = 1 for any Thus the 
lift-curve slope is 2ir for any Nc which is that predicted by thin wing theory. 


Chordwise Solutions for 

For Nc"^ and nc finite or Nc-Oc finite then equation (39) becomes 


Nc"“ 

n^=finite 

c 


'KK 
c c 


r(n^+l/2) 
,r(2n^-l) r(n^+T) 


Nc^ 

N,-n^=finite 
c c 




2(N^-nc+l/2) r(N^-n^+l/2) 
r(Nc-n V) ^ 


(48) 


For both and n^ large [using eq. (29)], equation (39) becomes 


N-xx. 


'N^n^ 
c c 


2(N -n 2 1-f 1/2 

= - (-^) 

\/C TT ' C ' 


TTfl, 


(49) 


Equation (49) is identical to the thin wing theory additional loading given 
inequation (40) and shows that a chord lattice solution converges to 
mathematical exactness as the lattice grid becomes infinite. With equations 
(48), (40), and (41) the chordwise loading factor at wing leading edge is 
given by (for Nc"^) 

f r(n-Ul) ^ ^ 

^ (50) 

c c (n^_i+i/4)i/^ r(n^-l+l/2) 

This function is identical to that of equation (31) when nc-1 = N-n, and is 
listed in equation (34). With equations (48), (40), and (41), the chordwise 
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loading factor at wing trailing edge is given by (for Nc^®) 


f (N -n +3/4)'/^ r(N -n tl) 

"c-Nc-V = IVn%r/^) WW^-n^ T /?) -- ' 

This function is identical to equation (33) when N^-n^ = N-n, and 

is listed in equation (34). 


Ag -VORTEX LATTICE 


Dependency of Spanwise Lattice on Effect of 
Aspect Ratio 

As aspect ratio approaches zero the spanwise optimum lattice arrangement 
is that defined by equation (3), that is at ''Ivn = (n-l/4)/N. As aspect 
ratio approaches infinity the spanwise lattice arrangement is the planform 
lattice which positions the trailing vortices at nvn = n/N, that is at the 
outside edge of the lattice panel including a vortex at the wingtip. This 
is because at A = ® the spanwise loading has the same distribution as the 
wing chord along the span. This high loading near the wing tip (when A = ®) 
must be taken into account by the lattice trailing vortices. The objective 
here is to develop an aspect ratio function factor for the lateral panel 
positions which asymptotically approaches the correct values at Ae-K) and at 
Ag-x”. The subscript e denotes the effective swept panel aspect ratio given 
by 


Ag = A/cosA-|/2 


(52) 


Tile planform lattice (nvn = i^/N) is the lateral lattice arrangement that 
has been in general use in most vortex-lattice methods. 


An aspect ratio equal to four is about an aerodynamic mean between 
0<Aei”. Loading values based on the Ag-^O lattice and on the Ae-x» lattice will 
be computed for an Ae = 4 rectangular wing. Comparison with the loading 
predicted by an accurate analysis will establish the A-effect function that is 
needed. By Biot-Sava'"t law the downwash due to an unskewed horseshoe vortex 
is 


W) 


nn. 


^0 ^A/To'-'^W^nXn^' ^0 V?o'''^’'wm“^An"\n^' 


m,m 


27iC, 


^wm'^An^'^vn 


(53) 


^wm'^An’^^vn 


Where = Cwmc " ?vnc» ^nnc the dimensionless vortex strength of the 
elemental horseshoe vortex. For Nc = 1 solution, Cq '"emains constant equal to 
1/4 for A = 4. lateral middle of the n panrl . For the slender 

wing (Ag->0) lattice use 
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= !!!:|a , "a, = 0 


For the planform lattice use 


= n = m-1/2 

’^vn ~ N * '^wm ~ N 




The equation to be solved for additional loadino for an Nq = 1 solution, 
A = 4 rectangular wing is 


N r 
n= 1 


Nn)^ V4VV16^(n^- \n^^ 


’^wm^'^vn 




Wing lift-curve slope is determined from equation (24) and wing loading 
by equation (20) with which 

^ i r 

^lSv “Nn 

where for the A-»0 lattice iw^ values are given in table 3, while for the 
planform lattice i|\|n are all unity. Results of the solutions of equation 
(56) with the Ae->0 lattice of (54) and with the planform lattice of (55), 
for increasing N, are presented in tables 6 and 7. An accurate loading 
prediction of this A = 4 rectangular wing is made in reference 19. From 
reference 19 

C, = 3.6623, n = cos^i 1 


= ^sind) + .07879sin3(t) + .01290sin5<j) +.00350sin7({) + .00170sin9(|)) J 

IT 

Percent difference*" from the values of equation (58) are shown in tables 6 and 
7. For this A = 4 wing the A->0 lattice has less than one-half the error of 
the planform lattice. 

Spanwise A^-Vortex Lattice 

Since for the A = 4 wing the A-K) lattice prediction is too small and the 
planform (A-ho) lattice prediction too large then a factor can be developed 
between the two w hich forms the bas-is of spanwise Ae-lattice arrangement. 

This factor is 2/VAe+4 , then for the Ae-lattice 
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vn 


'ANn 


_ 1/ 1 X 1. 1 1 

^ 24 K~H ’ ^ N I 




+4 
e 


>(59) 


1^1 1 ^ - 1 In 

m(tt - -==^ . \ " n(n- J =r) , 

n iM ^ J 


where Ap is defined in equation (52) and is the spanwise position of the 
computed loading distribution. These spanwise lattice panel positions 
asymptotically approach the A->0 lattice and the planform lattice as A^D and ~ 
respectively. With equation (59) the solutions of equation (56) for lift and 
loading are listed in tables 6 and 7. With the Ae-lattice the predicted lift 
and loading is accurate for all N's which was the basis for the term 
optimum applied to the A->0 lattice and chordwise lattice. Equation (59) is 
simply a mathematical statement tha t re lative to the planform lattice all the 
panels are shifted inboard by l/2NVAe+4, and that the downwash point is at the 
lateral center of each panel except in the wing root panel. 


Application of the A^-Vortex Lattice 

For a unifojT vortex lattice, the elemental skewed horseshoe vortex 
lateral position at the trailing vortices is given by nvn in equation (59) and 
the chordwise position by Cvnc equation (36). The downwash points are 
positioned laterally at nwm in equation (59) and chordwise at Cwmc in 
equation (36). Let Tnnc unknown circulation strength of the elemental 

skewed horseshoe vortex. For symmetric wing loading rnnr is determined from 
an NNc matrix solution which satisfiesNNc downwash point boundary conditions. 
The pressure coefficient at span station nn in equation (59) and chord 
station = (nc~3/4)/Nc is given by 


= -2N i f 

i^nn ^^c^ANn'N n.-— rr 
c c c c V 

n 


(60) 


where spanwise and chordwise loading factors are included [eq. (59), tables 3 
and 5]. Now 


r 


n 



then the spanwise loading at of equation (59) is 


(61) 


/ 

'ias- ,, i 
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atn ANn n 
n 


(62) 


\/hich includes the spanwise loading factor [eq. (59), table 3]. The spanwise 
loading gradient at span station nvn equation (59) is 


dn" ' ‘^ANn ^^^n'Vl^ 


(63) 


where rn is given in equation (61) and where the spanwise loading gradient 
factor is given by 


'’ANn = '•''■"Nn'^ 


, with h|^^ in table 2 . 


The lift coefficient is [r^ from eq. (61)] 

'^L ' ^1''' * Z ''n> 


N 

z 

n=2 


(64) 


The induced drag is given by 

C = -- 
''D. 2bV 


1 N 

/ “i'''’’' ‘ * Z Vn> 

' n-2 


- JL 

'n " TrbV 


N-1 

1 1 ** 1 

^N-n+l/2'1r+?rT^^N z, ^m^+l/2Vn+l^^ WlM 

*“ m=l 


(65) 


where Tp is given in equation (61). For spanwise loading distribution already 
known then equation (65) provides a convenient method for predicting induced 
drag. For this case rp)gtnp is known then Fp is determined from equation (62) 
for application in (65). 


CONCLUDING REMARKS 


The application of the vortex-lattice method to the slender wing configu- 
ration has provided a rigorous analytical basis for the spanwise properties 
of the vortex-lattice method. Mathematical similarities are shown between 
the spanwise panel lattice solution and the chordwise solution which con- 
verge to thin winq theory results. As the number of chordwise panels approaches 
infinity, thin wing theory chordwise loading is predicted exactly except in 
the limit points exactly at the leading edge and at the trailing edge. At 
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these two points chordwise loading factors are mathematically evaluated which 
are useful in finite panel solutions. Similarly, as the number of spanwise 
panels becomes infinite, slender wing theory spanwise loading is predicted 
exactly except at the point of the wing tip. At this point a spanwise loading 
factor is mathematically determined from a limit solution. The presentation 
in this paper is based on a planform uniform distribution of panels chordwise 
and spanwise. In a discussion in 1972 Mr. W. B. Kemp, Jr. of NASA-Langley 
said he had found that the chordwise vortex-lattice solution gave an accurate 
integrated lift for an arbitrary planform panel distribution along the chord. 
As part of the present work this was investigated and it was shown that using 
the 1/4 - 3/4 rule for locating the vortex and downwash point in the planform 
panel, the chordwise lift-curve slope of 2tt and also the spanwise slender wing 
value of ■nk/2 are predicted for any distribution of planform panels on the 
wing and for any total number of panels. However, the loading distribution 
factors are not as near unity. The aspect ratio effect on spanwise lattice 
arrangement has a weak chordwise counterpart for a chordwise lattice 
arrangement. An initial study of this effect indicates that the chordwise 
lattice arrangement differs when aspect ratio is less than unity. In 
conclusion, the Ag-vortex-lattice arrangement of the previous section provides 
a uniform uncomplicated, accurate system. It leads to computations with a 
high accuracy to work ratio. 
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TABLE 1. - SPANWISE SOLUTIONS OF EQUATION (6) 


\n 
n\ 1 

2 

3 

4 

5 6 

1 1 
2 M 

3x5 

2x3*+ 

3x5x7 




3 33 

23x3** 

2x3**x53 



3x5x7 

3x5x7x9 

llodd! 



a. 22x32x5 
3x5x7x9 

22x33x5 

22x33x53 
1 3odd • 

22x33x53x73 
1 5odd • 


i; 32x52x7 

23x3^x52 

l3odd* 

2x33x53x7 

'5odd! 

23x33x53x73 
^ ^odd' 

2x33x53x73x93 

l^odd! 

c 2x3^x5x7 

® TW 

1 

2x3^x52x72 
^ ^odd • 

2S35x5“x7 

2**x33x53x73 
^ ^odd • 

22x33x5**x73 22x33x52x73x 11 3 
^^odd* \^^dd- 

3x 2N-1 [(2N-l)odd!]^ 
(N-l)!(4N-1)odd! 

T 

3N'^N-3)odd! 3N 




TABLE 2. - SPANWISE LOADING GRADIENT FACTOR 



1 

2 

3 

4 

5 

6 

1 

1 .082788 






2 

.965721 

1 .118660 





3 

.958851 

1.000148 

1 .124345 




4 

.956981 

.993739 

1 .005515 

1.126193 



5 

.956199 

.992082 

.999207 

1.007242 

1 .127013 


6 

.955795 

.991404 

.997613 

1.000965 

1 .008002 

1.127446 

oo ^ 

.954930 

1 

.990297 

.995956 

.997804 

\ 

.998625 

\ 

.999060 

-N+r. 
00 ^ 

1 .000564 

1 .000858 

1 .001459 

1.003004 

\ 

1 .009253 

\ 

' . 1 28379 
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TABLE 3. - 
2 

SPANWISE LOADING FACTOR 
3 4 5 

6 

N 

11=0 






1 

.954930 






2 

.983016 

.966313 





3 

.990965 

.988798 

.972705 




4 

.994332 

.993811 

.992034 

.974770 



5 

.996085 

.995935 

.995712 

.993301 

.975684 


6 

.997119 

.997082 

.997167 

.996548 

.993923 

.976166 


-N+n 
+6 


.999530 


.999312 .998901 


.997976 


.995137 


.977205 


TABLE 4. - e , CHORDWISE SOLUTIONS OF EQUATION (3>) 
t c 

K K 

c c 


n 1 

NX 

2 

3 

4 

5 

® n -l\"r 

n^-1 c c 
c 

/ . (n 

n =1 " ' 

c 

C 

1 1 





1 

1/4 

2 3/2 

1/2 




2 

4/4 

3 15/8 

3/4 

3/8 



3 

9/4 

4 35/16 

15/16 

9/16 

5/16 


4 

16/4 

5 315/128 

35/32 

45/64 

15/32 

35/128 

5 

25/4 

6 693/256 

315/256 

105/128 

75/128 

105/256 

63/256 6 

36/4 






II 

II 






N 

N-/4 


)e^ 


c c 
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1.102658 
1 .122892 


.986247 


1.126095 1.004345 
1.127139 1.007210 


.004345 .980140 

.007210 .998126 

.008144 1.000974 


1.127603 1.008144 1.000974 
1.127848 1.008559 1.001902 


.978631 

.996589 

.999433 


.978043 

.995991 .977756 


\ \ \ \ 

1.128379 1.009253 1.003004 1.001459 1.000858 1.000564 


.999530 .999312 .998901 .997976 .995137 .977205 



TABLE 7. - LATTICE COMPARISONS FOR PREDICTING X A= 4 RECTANGULAR WING 

^L 



planform lattice 


A 

-<0 lattice 



A ■ 
e 

•lattice 




H 






Cl 




c 

1 

. Jw n/ 

i 

N 

n 


^L,19 

Ap— ,% 

n 


<=L 

'T.is 

Ap »% 

^L 


^L 

f 

L,19 

As— 

^L 

ANn 

1 

.500 

1 

1 .092 

-8.45 

0 

T 

.273 

1.187 

7.25 

.146 

1 .176 

T 

.179 

-. 29 

.9681 

2 

.250 

1 .088 

1.165 

-6.62 

0 

1 

.230 

1.187 

3.65 

.073 

1 .185 

1 

.185 

0 

.9880 


.750 

.912 

.914 

- .21 

.625 

1 

.026 

1.024 

.19 

.662 

.989 


.995 

- .67 

.9762 

3 

.167 

1 .126 

1.177 

-4.33 

0 

1 

.224 

1 .187 

3.09 

.049 

1 .194 

1 

.187 

.58 

.9936 


.500 

1.049 

1 .092 

-4.01 

.417 

1 

.145 

1 .124 

1 .87 

.441 

1.115 

1 

.115 

- .04 

.9921 


.833 

.826 

.799 

3.34 

.750 


.891 

.914 

-2.57 

.77’ 

.871 


.883 

-1.36 

.9807 
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Figure 1.- Comparison of chordwise loading with spanwise loading gradient 




